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ABSTRACT 



Context. NRAO 150 -a compact and bright radio to mm source showing core/jet structure- has been recently identified as a quasar at 
redshift z = 1.52 through a near-IR spectral observation. 

Aims. To study the jet kinematics on the smallest accessible scales and to compute the first estimates of its basic physical properties. 
Methods, we have analysed the ultra-high-resolution images from a new monitoring program at 86 GHz and 43 GHz with the GMVA 
and the VLBA, respectively. An additional archival and calibration VLBA data set, covering from 1997 to 2007, has been used. 
Results. Our data shows an extreme projected counter-clock-wise jet position angle swing at an angular rate of up to « 1 l°/yr within 
the inner « 31 pc of the jet, which is associated with a non-ballistic superluminal motion of the jet within this region. 
Conclusions. The results suggest that the magnetic field could play an important role in the dynamics of the jet in NRAO 150, which 
is supported by the large values of the magnetic field strength obtained from our first estimates. The extreme characteristics of the jet 
swing make NRAO 150 a prime source to study the jet wobbling phenomenon. 
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1. Introduction 

An increasing number of jets in active galactic nuclei (AGN) 
have been reported to show either regular or irregular swings 
of the innermost jet structural position angle in the plane of the 
sky (e.g., in OJ 287, Tateyama & Kingham 2004; in 3C 273, 
Savolainen et al. 2006; in 3C 345, Lobanov & Roland 2005; in 
BL Lac, StirUng et al. 2003; in S5 0716+71, Bach et al. 2005). 
Time scales between 2 and 15 years and structural position-angle 
oscillations with amplitudes from ~ 25° to ~ 45° are typical for 
the reported cases. We will call this phenomenon jet wobbling 
hereafter. Parsec scale AGN jet curvatures and helical-like struc- 
tures also at larger distances from the central engine are also be- 
lieved to be triggered by changes in direction at the jet ejection 
nozzle (e.g., in 3C 84, Dhawan et al. 1998). 

The physical origin for the observed jet wobbling is still 
poorly understood. Among the various possibilities, regular pre- 
cession of the accretion disk is frequently used for modeling 
at present. Most AGN precession models are driven either by 
a companion super-massive black hole or another massive ob- 
ject inducing torques in the accretion disk of the primary (e.g.. 
Lister et al. 2003, for 4C +12.50; Stirling et al. 2003, for BL Lac; 
Caproni & Abraham 2004 for 3C 120) or by the Bardeen- 
Peterson effect (e.g., Liu & Melia 2002; Caproni et al. 2004). 
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However, other AGN scenarios that have yet to be explored ex- 
tensively, such as the orbital motion of the jet nozzles (also in- 
volving binary systems) or other kinds of more erratic disk/jet 
instabilities (e.g., similar to those thought to produce the quasi 
periodic oscillations [QPO] in X-ray binaries), can not be ruled 
out yet. Note that, in support of these erratic instabilities, it is 
still under debate whether the observed jet wobbling is strictly 
periodic or not (see Mutel & Denn 2005 for the case of BL Lac). 

There is still no paradigm to explain the phenomenon of jet 
wobbling in AGN, but it is rather likely that, as it is triggered in 
the innermost regions of the jets, it must be tied to fundamental 
properties of the inner regions of the accretion system. Hence, 
there is ample motivation to study the jet wobbling phenomenon 
to place our understanding of the jet triggering region and the 
super-massive accretion systems on firmer ground. 

VLBI observations at millimetre wavelengths are a powerful 
technique to image the innermost regions of AGN jets -which are 
self-absorbed at longer wavelengths- with the highest angular 
resolutions; ~ 50//as at 86 GHz (3.5 mm) and ~ 0.15 mas at 
43 GHz (7 mm). Here, we report the discovery of an extreme 
case of jet swing in the quasar NRAO 150, through the first ultra- 
high-resolution VLBI set of images obtained from this source at 
86 GHz and 43 GHz. NRAO 150 is a strong radio-mm source. 
At radio wavelengths, on VLBI scales, NRAO 150 displays a 
compact core plus a one-sided jet extending up to r i 80 mas 
with a jet structural position angle (PA) of ~ 30° (e.g., Fey & 
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Chariot 2000). NRAO 150 was not detected by the early optical 
surveys, most probably due to obscuration through the Milky 
Way (Galactic latitude = -1 .6°). Almost 40 yr after its discovery 
at radio wavelengths (Pauliny-Toth, Wade & Heeschen 1966), 
NRAO 150 has been identified as a quasar at redshift z - 1.52 
through a near-IR spectroscopic project (Acosta-Pulido et al. in 
prep.). 

Throughout this paper we assume Ho - 72 km s"' Mpc 
Q.m = 0.3, and £1a = 0.7. Under these assumptions, the lumi- 
nosity distance of NRAO 150 is di = 11025 Mpc, 1 mas corre- 
sponds to 8.5 pc in the frame of the source, and an angular proper 
motion of 1 mas/yr translates into a speed of 69.4 c. These pa- 
rameters are used here for the first time to make quantitative es- 
timates of the basic physical properties of the jet in NRAO 150. 



2. Observations, images and their modelling 

The data set presented here consists of 5 VLBI images at 
86 GHz taken with the Global Millimetre VLBI Array (GMVA', 
and its predecessor, the CMVA) covering the time range from 
October 2001 to April 2004, and of 34 Very Long BaseHne Array 
(VLB A) images taken at 43 GHz from June 1997 to February 
2007. Among these 34 images, 20 of them (those obtained from 
June 1997 to May 2002) were taken as calibration measurements 
of the SiO maser monitoring program of TX Cam (Diamond & 
Kemball 2003). A new 8.4 GHz VLBA image obtained in June 
2002 is also presented here. Information about each single ob- 
servation and its resulting image is provided in Table 1 . 

The 86 GHz and the 43 GHz data were correlated at the 
Max-Planck-Institut fiir Radioastronomie and the VLBA corre- 
lators, respectively. Subsequent phase and amplitude (including 
opacity) calibration were performed using standard procedures 
within the AIPS software. Further phase and amplitude self- 
calibration, as well as the final imaging, were performed through 
standard procedures within the Difmap software. 

Examples of the resulting images at 86 GHz, 43 GHz, and 
8.4 GHz -obtained in mid 2002- are presented in Fig. 1 to- 
gether with a 2.3 GHz VLBA image obtained in July 2005 by 
Kovalev et al. (2007). At the four observing frequencies, the im- 
ages show the typical core plus one sided jet structures of radio- 
loud quasars. However, the orientation of the sub-mas scale jet 
(Fig. 1-a and b) diff'ers by > 100° with respect to that at larger 
scales (Fig. 1-c and d). This suggests a bent structure of the inner 
jet oriented within a very small angle to the line of sight. 

The series of 43 GHz images in Figure 2 shows evidence 
of motion in the North- West to South-East direction. More in- 
triguing is the fact that this motion is accompanied by an evident 
counter-clock-wise rotation of the jet on the plane of the sky, 
which will be the subject of discussion hereafter 

Difmap was also used to model the images of NRAO 150 as 
sets of 3 or 4 circular-Gaussian intensity distributions (model 
components, e.g.. Fig. 1). Elliptical Gaussian fits were also 
tested, but they did not improve the results. Only the minimum 
number of model components describing jet regions consider- 
ably brighter than the image-background rms and at close posi- 
tions during contiguous epochs were allowed. All fitted model 
components resulted to have flux densities 5'comp(86 GHz) > 
90mJy, 5comp(43GHz) > 50mJy, 5comp(8.4GHz) > 60mJy. 
This characterized reliably each of the most significant jet re- 
gions by its flux density, relative position, and a FWHM size 
measurement at each observing epoch. Table 2 lists the results 
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Fig. 1. VLBI images of NRAO 150 at 86 GHz (a), 43 GHz (b), 8.4 GHz 
(c), and 2.3 GHz (d). The 2.3 GHz image was obtained by Kovalev 
et al. (2007). In images a, b, c and d, the lower contours correspond 
to 1.0%, 0.2%, 0.15% and 0.5% of the corresponding image peak of 
1.05 Jy/beam, 3.84 Jy/beam, 5.55 Jy/beam, and 0.69 Jy/beam, respec- 
tively. The upper contours increase by factors of 2. Other image param- 
eters of the a), b) and c) maps are given in Table 1. The 2.3 GHz con- 
volving beam is (7.4 x 2.6) mas^ with major-axis position angle at -2°. 
The 86 GHz image illustrates the capability of the GMVA to achieve 
angular resolutions of up to 40 yuas with dynamic ranges of up to 100:1 
and notable image fidelity -as suggested by the matching 86 GHz and 
43 GHz source structures-. In the 43 GHz image, the circles are centred 
at the positions of the Gaussian model components. The radius of such 
circles symbolize the FWHM size of the corresponding model compo- 
nent. 

from such fits, using the position of the (typically) brightest 
modelled Gaussian component -the VLBI core- as reference. 

3. Superluminal motion and non-ballistic jet swing 

The 43 GHz jet structure evolution is illustrated in Fig. 2, 
whereas Fig. 3 shows the 86 GHz and 43 GHz inner jet axis 
evolution, the projected trajectories of the three identified model 
components Ql, Q2, and Q3 with respect to the core position, 
and the total flux evolution of these components. We modelled 
these trajectories, taking into account their curvatures, by fit- 
ting second order polynomials to them (e.g., as in Homan et 
al. 2001 and Jorstad et al. 2005); see the results in Table 3. 
The mean measured proper motions are (0.047 + 0.002) mas/yr, 
(0.042 + 0.001) mas/yr, and (0.033 + 0.002) mas/yr, for Ql, Q2 
and Q3, with corresponding superluminal apparent speeds (J3°^p) 
of(3.26+0.14)c, (2.85+0.07)c, and (2.29+0. 14) c, respectively. 
Fig. 3 shows that the motion of features Ql, Q2 and Q3 outwards 
from the core is accompanied by counter-clock-wise rotation of 
their position vector, projected onto the plane of the sky, with 
extremely fast angular speeds < 0(Q1) >= (10.7 + 0.7)°/yr, 
< 0(Q2) >= (5.8 + 0.3)7yr and < 0(Q3) >= (8.0 + 0.5)7yr. 

These are the first superluminal proper motions reported 
within the innermost 0.5 mas of the jet in NRAO 150. We 
should also emphasize that the three fitted trajectories are bent 
and non radial, which implies that the motion of all the jet fea- 
tures are non-ballistic up to a projected distance from the core of 
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~ 0.5 mas. Assuming a jet viewing angle <p ^ 8° (see Section 4), 
the jet features are implied to behave non-balhstically up to a 
deprojected distance from the core « 31 pc. 

Similar results, in terms of the magnitude of projected angu- 
lar speeds and non-ballistic nature, are obtained if either Ql or 
Q2 is selected as the kinematic centre of the source. 

Note also that when the mean projected velocity vectors of 
Ql, Q2 and Q3 are decomposed in their radial (with respect to 
the core) and non-radial (i.e. normal, y6°pp_norm) components, the 
later are (2.71 + 0.21) c, (1.39 + 0.21) c, and (2.01 + 0.14) c, re- 
spectively. Thus, we have shown here an extreme case of super- 
luminal non-ballistic swing in the jet of an AGN. This is unusual 
because superluminal motions are reported in most cases in the 
radial direction relative to the core only. Recently however, a 
growing number of AGN jets exhibiting bent model-component 
trajectories have been reported (e.g., Homan et al. 2001; Jorstad 
et al. 2005). NRAO 150 seems to be an extreme case, which 
reflects the remarkable non-ballistic nature of its jet. Both this 
non-ballistic nature and projection effects can explain why the 
trajectories of Ql and Q3 do not cross with the one of Q2 when 
they are extrapolated back (see Fig. 3). 

4. Jet physical parameter estimates 

Figure 3 shows the total flux density evolution of the model com- 
ponents. Dramatic increases are displayed by Q2, at ~ 2000 and 
~ 2005, and by Q3, at ~ 2002 and ~ 2006. Interestingly, a pro- 
nounced dip in the core light curve happens shortly before each 
of the dramatic Q2 flux peaks. These changes can be explained 
by changing Doppler boosting factors due to either a jet feature 
propagating through a bent underlying jet, or to the propagation 
of such a bend itself along the jet. 

The "synchrotron-self-Compton" (SSC) Doppler factor of 
the source ((5ssc) can be constrained by comparing the observed 
and the expected first order self-Compton X-ray flux density 
with the flux density of the source (Sv,„) at the synchrotron 
turnover frequency (v„,; e.g., Marscher 1983; Agudo et al. 2006). 
Between August 1990 and February 1991, NRAO 150 displayed 
flux densities of Sikev = (0.35 + 0.09) /jJy at 1 keV. ROSAT 
measured a flux /^o. 1-2.4 keV = (0.22 + 0.06) x 10"" erg cm"^ s~^ 
corrected for Galactic absorption (assumed a photon index F = 
-1.7 + 0.1) from IRXS J035930.6-H505730 (Voges et al. 2000), 
which position is fuUy consistent with the one of NRAO 150 
even at the level of the ROSAT l-cr position error. No other 
X-ray source or bright radio-mm source is found within 1 ar- 
cminute of the position of NRAO 150. At * 37 GHz, the flux 
density Sy^ = (2.59 ± 0.21) Jy (see Reuter et al. 1997; Voges et 
al. 2000; Chen et al. 2001; Terasranta et al. 2005). Given that flie 
period between 1990 and 1997 was characterized by a low ac- 
tivity phase and that the VLBI structure of the source was dom- 
inated by the core in the 1990s, we assume that the bulk of the 
37 GHz flux measured at the beginning of the 1990s was radiated 
from the core. By taking into account these measurements, the 
average 43 GHz size of the core in our VLBI maps from 1997 
to 2007, < e >= (0.10 ± 0.03) mas, and the typical optically 
thin synchrotron spectral index a » -0.5 at the beginning of the 
1990s, as well as the assumptions and approximations outUned 
in Agudo et al. (2006), we obtain 5ssc ~ 6. 

By using expression (2) in Marscher (1983) and the outlined 
measurements for the core (5v,„ ~ 2.59 Jy, ^ 37 GHz, 6 x 
0.10 mas, ^ssc ~ 6 and a ^ -0.5), the magnetic field intensity 
in the jet can be estimated. In this way we find B « 0.7 G, which 
is rather large for an AGN jet with the X-ray emission produced 
by SSC (e.g., Marscher et al. 1979; Ghisellini et al. 1998). 



Using now both /3^^(Ql) as an estimate of the jet plasma 
speed and 6 » die in the definitions of and the apparent proper 
motion of the fluid, /3app^, we obtain y a; 4 and a: 8°. 

Note that all these estimates should still be taken as first 
approximations, since the used measurements are not as accu- 
rate and simultaneous as needed. The results from a coordinated 
multi-waveband campaign organized in February 2007 will al- 
low better and updated estimates. In addition, the nature of Ql, 
Q2 and Q3 is still not known. They might be related to pattern 
speeds of propagating Kelvin-Helmholtz perturbations, which 
must have smaller speeds than those of the fluid in which they 
propagate (e.g., Perucho et al. 2006). In contrast, they could be 
shock waves propagating in a bent underlying jet. In this case, 
the speed of the shocks must be an upper hmit of the plasma 
speed. Nevertheless, our estimates of y and (p are consistent 
with those expected from the iimermost regions of a quasar like 
NRAO 150, which extreme (projected) jet misaUgnmentimpUes 
a very small viewing angle. 

5. Summary and conclusions 

We have reported the results from the first 86 GHz and 43 GHz 
VLBI monitoring program of the recently identified quasar 
NRAO 150. Our observations reveal a) a large projected mis- 
alignment of the jet by > 100° within the inner 0.5 mas to 1 mas 
from the core, b)an extremely fast counter-clockwise rotation of 
the projected jet axis at a rate of ~ 6°/yr to ~ ll°/yr, c) non- 
ballistic superluminal motions with mean speeds from 2.3 c to 
3.3 c within the inner 0.5 mas from the core (deprojected dis- 
tance X! 31 pc), d) transverse (non-radial) speeds of 2.7 c, 1.4 c, 
and 2.0 c for Ql, Q2 and Q3, respectively, e) an extreme case of 
superluminal non-ballistic jet swing, and/j the first approxima- 
tions to quantitative estimates of the basic physical properties of 
the jet in NRAO 150: 6 ^ 6, B ^ 0.7 G, y ^ 4, and (p * 8°. 

Whereas the ultimate origin of the jet swing must be an in- 
trinsic change of the direction of the iimer jet axis (either caused 
by changes produced at the injection region or by interaction 
with the medium surrounding the jet), possible causes for the 
non-ballistic nature of the emitting flow are either a small iner- 
tia of the jet compared to that of the impacted ambient medium, 
or by the fact that we are observing a jet instability propagating 
downstream. However, the superluminal non-ballistic motion of 
the jet features in NRAO 150 is perhaps too fast and system- 
atic along tens of parsecs (deprojected, see above) to be induced 
solely by the ambient medium. In addition, the perturbation that 
must be produced either by the impact with such medium or by 
the changes at the injection region should imply the growth of 
disruptive instabilities (e.g., Perucho et al. 2006), in contrast to 
the remarkable colhmation of the jet up to the kiloparsec scale 
(see Fig. 1). Mizuno et al. (2007) have shown that a magnetic 
field with the appropriate configuration helps keeping a jet colli- 
mated against the growth of instabilities. Hence, this possibility, 
together with the high value for the magnetic field intensity es- 
timated in Section 4, suggests that the magnetic field could play 
an important role in the dynamics of the jet in NRAO 150. 

It is still unclear whether the reported change of the direction 
of ejection in NRAO 150 is related to a regular (strictly periodic 
or not) behaviour or to a single event. This, together with the still 
unknown nature of the moving jet features (either propagating 

^5 = [y(l - ficos(p)r^ where y = {1 - jf)'^'^ is the Lorentz factor, p 
is the speed in units of the speed of light and 4> is the angle between the 
direction of the flow and the line of sight 

^ jSapp = (fisin<p)/{l -ficos<p) 
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Fig. 2. Selection of 18 of the 34 NRAO 150 images obtained from 1997 to 2007 with the VLB A at 43 GHz. All images have been convolved 
with the same circular beam of 0.16 mas. The common intensity scale ranges from 0.025 Jy/beam to 4,128 Jy/beam. The positions of the fitted 
Gaussians are indicated by the black crosses, whereas the black circles (of radius equal to the FWHM of each Gaussian) symbolize their size. A 
movie corresponding to these images can be downloaded from <http://www.iaa.es/~iagudo/research/NRAO150/NRAO150VLBA43GHznc.avi>. 



86GHz/45GHz jet axis evoljtion 43 GHz poslTions evolution 43 GHz flux evolution 




Rel.' R.A. (mos) ' ' ' Rel. R.A. (mos) '938 2000 2002 2004 2006 2008 1998 2000 2002^ ^2004 2006 2008 

Fig. 3. From right to left, 86 GHz and 43 GHz inner jet axis evolution, projected trajectories of the 43 GHz model components, and their cor- 
responding 43 GHz flux evolution. The crosses in the middle plot symbolize the position errors (computed using the prescriptions by Jorstad et 
al. 2005), whereas the superimposed bent lines represent the fitted trajectories of each model component (labelled in the Figure). Filled diamonds 
in the flux evolution plots symbolize upper limits (see Table 2). 



curvatures or shocks in a bent jet), does not allow us to specify 
the ultimate origin of this phenomenon in the source. This will 
be the matter of future studies. Nevertheless, the extreme char- 
acteristics of the jet swing make NRAO 150 a prime source for 
future studies of the origin of jet wobbling. 
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Table 1. 86 GHz and the 43 GHz image information. Given are the observing epochs, the total integration times (fi„t), the observing frequency 
bandwidth (Avobs), the number of bits used for the signal digitalization sampling, the FWHM minor and mayor axes of the restoring beam and its 
orientation angle, the integrated total flux densities, the peak flux density, and the noise levels of the resulting images. 
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20 


16 


2 


129 


241 


-19.0 


5.73 


3.30 


4.0 


2000.72 


20 


16 


2 


132 


244 


-19.0 


6.04 


3.39 


3.5 


2000.80 


20 


16 


2 


132 


266 


-21.3 


6.18 


3.42 


3.4 


2000.88 


20 


16 


2 


242 


356 


-36.7 


6.18 


4.75 


4.4 


2000.97 


20 


16 


2 


205 


329 


23.4 


5.92 


3.85 


4.5 


2001.05 


20 


16 


2 


208 


340 


13.8 


5.72 


3.85 


4.1 


2001.13 


20 


16 


2 


206 


268 


5.7 


5.74 


3.40 


3.7 


2001.20 


20 


16 


2 


195 


261 


8.8 


5.73 


3.35 


5.8 


2001.29 


20 


16 


2 


225 


275 


12.7 


5.59 


3.45 


4.2 


2001.38 


20 


16 


2 


1 j4 




-19.5 


5.40 


2.49 


3.1 


2001.48 


20 


16 


2 


138 


238 


-23.6 


5.31 


2.57 


3.2 


2001.62 


20 


16 


2 


137 


256 


-31.7 


5.41 


2.83 


4.5 


2001.70 


20 


16 


2 


129 


239 


-28.5 


4.79 


2.28 


3.5 


2001.78 


20 


16 


2 


134 


248 


-21.3 


5.69 


2.63 


3.6 


2001.86 


20 


16 


2 


135 


243 


-20.7 


6.76 


3.58 


2.7 


2002.07 


20 


16 


2 


160 


249 


-20.8 


7.74 


3.99 


2.9 


2002.40 


25 


16 


2 


141 


212 


-16.0 


7.75 


4.27 


3.0 


2002.48 


220 


32 


2 


139 


184 


-3.6 


7.60 


3.92 


0.8 


2003.20 


220 


32 


2 


147 


186 


1.1 


6.19 


2.60 


0.9 


2003.49 


220 


32 


2 


140 


186 


0.1 


7.36 


3.42 


0.8 


2004.64 


220 


32 


2 


119 


170 


-8.3 


7.02 


3.22 


1.1 


2004.97 


220 


32 


2 


129 


196 


-26.5 


5.37 


2.10 


0.9 


2005.25 


220 


32 


2 


121 


177 


3.1 


4.82 


1.59 


1.0 


2006.17 


84 


64 




114 


190 


2.7 


4.27 


1.91 


1.5 


2006.24 


84 


64 




113 


207 


1.4 


4.31 


2.01 


1.5 


2006.30 


84 


64 




111 


186 


-1.0 


4.34 


1.91 


1.6 


2006.37 


84 


64 




110 


187 


-3.8 


4.20 


1.77 


1.9 


2006.46 


84 


64 




112 


188 


-14.8 


4.21 


1.83 


1.9 


2006.88 


42 


32 


2 


130 


276 


-32.4 


4.26 


2.33 


1.2 


2007.16 


42 


32 


2 


124 


266 


-31.4 


6.32 


3.07 


1.7 


8.4 GHz images 


2002.48 


84 


32 


2 


1180 


904 


-3.6 


6.73 


5.65 


0.9 
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Table 2. 86 GHz, 43 GHz and 8.4 GHz circular-Gaussian model-fit parameters. S ,r,@ and FWHM are the flux density, the projected distance to 
the Core, the position angle respect to the Core and the FWHM size, respectively. Upper limits on the FWHM size and S corresponds to fits of 
unresolved jet regions. For such fits, the upper limit of the sizes were imposed in Difmap to estimate the corresponding flux density upper limits. 
The latter are close approximations to the actual flux densities of the unresolved regions. 



Comp. 


s 


r & 


FWHM 




(Jy) 


(mas) n 


(mas) 


86 GHz -2001.82 


Core 


1.937+0.019 




0.053+ 0.001 


WO 


0.661+0.020 


0.13+0.02 -109.9+9.0 


0.137+0.007 


W2 


< 0.244 


0.28+0.01 -140.5+2.1 


<0.010 


W3 


< 0.294 


0.22+0.01 -130.8+2.6 


<0.010 


86 GHz - 2002.30 


Core 


1.799± 0.018 




0.042+ 0.000 


WO 


0.470+ 0.047 


0.18+0.02 -129.2+5.3 


0.034+ 0.003 


W2 


< 0.111 


0.31+0.01 -146.2+1.8 


<0.010 


W3 


< 0.159 


0.26+0.01 -139.2+2.2 


<0.010 


86 GHz -2002.81 


Core 


2.332± 0.023 




0.050+ 0.001 


WO 


0.933± 0.028 


0.15+0.02 -132.9+7.8 


0.106+ 0.005 


W2 


0.278+ 0.028 


0.32+0.01 -149.0+2.1 


0.024+ 0.002 


W3 


< 0.138 


0.22+0.01 -145.8+2.6 


<0.010 


86 GHz -2003.32 


Core 


1.376+ 0.014 




0.059+ 0.001 


W2 


< 0.224 


0.30+0.01 -151.6+1.9 


<0.010 


W3 


0.326+ 0.033 


0.22+0.02 -146.1+4.8 


0.037+ 0.004 


86 GHz - 2004.29 


Core 


3.663+ 0.037 




0.092+ 0.001 


W2 


0.093+ 0.009 


0.37+0.01 -152.4+ 1.7 


0.021 + 0.002 


W3 


0.785+0.024 


0.25+0.02 -147.0+4.6 


0. 111 + 0.006 


43 GHz- 1997.51 


C 


3.046+ 0.253 




0.085+ 0.001 


Qi 


< 0.180 


0.18+0.03 22.2+ 10.2 


<0.030 


Q2 


< 0.061 


0.18+0.03 132.6+10.3 


<0.030 


43 GHz -1998.14 


C 


3.051+0.147 




0.073+ 0.001 


Ql 


0.573+ 0.032 


0.13+0.03 43.4+ 13.8 


0.111+0.006 


Q2 


< 0.485 


0.12+0.01 118.7+4.9 


<0.030 


43 GHz - 1998.88 


C 


3.855+0.186 




0.106+0.005 


Ql 


i 0.154 


0.17+0.03 32.5+ 10.5 


<0.030 


Q2 


< 0.224 


0.14+0.03 122.8+ 12.9 


<0.030 


43 GHz - 1999.10 


C 


3.534+ 0.120 




0.066+ 0.001 


Ql 


0.602+ 0.027 


0.13+0.03 50.5+ 14.0 


0.101+0.005 


Q2 


< 0.653 


0.13+0.01 121.1+4.4 


<0.030 


43 GHz - 2000.54 


C 


3.550+ 0.243 




0.084+ 0.001 


Ql 


0.685+0.051 


0.13+0.03 49.3+ 13.8 


0.136+ 0.007 


Q2 


1.480+ 0.101 


0.15+0.01 128.5+3.7 


0.057+ 0.001 


43 GHz - 2000.72 


C 


3.645+0.154 




0.087+0.001 


Ql 


0.619+0.031 


0.15+0.03 50.2+ 12.5 


0.133+0.007 


Q2 


1.794+ 0.076 


0.16+0.01 130.5+3.5 


0.078+ 0.001 


43 GHz - 2000.80 


C 


3.562+ 0.187 




0.085+ 0.001 


Ql 


0.838+ 0.050 


0.13+0.03 53.7+ 14.1 


0.160+ 0.008 


Q2 


1.744+ 0.091 


0.17+0.01 130.4+3.3 


0.080+ 0.001 
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Table 2. continued. 



Comp. 


S 


r 





FWHM 




(Jy) 


(mas) 


n 


(mas) 


43 GHz - 2000.88 


C 


3.801+0.177 






0.088+ 0.001 


Qi 


< 0.430 


0.17+ 0.03 


60.1+ 11.1 


<0.030 


Q2 


1.944+ 0.106 


0.18+ 0.03 


133.6+ 10.1 


0.101+0.005 


43 GHz - 2000.97 


C 


3.885+0.145 






0.092+ 0.001 


Ql 


< 0.576 


0.16+ 0.03 


70.2+ 11.6 


<0.030 


Q2 


1.461 ±0.055 


0.20+ 0.01 


138.4+ 2.9 


0.067+ 0.001 


43 GHz -2001.05 


C 


3.618+0.111 






0.085+0.001 


Ql 


0.541+0.023 


0.17+ 0.03 


56.4+ 11.1 


0.058+ 0.003 


Q2 


1.581+0.048 


0.18+ 0.01 


131.6+3.1 


0.055+ 0.001 


43 GHz -2001.13 


C 


2.889+ 0.085 






0.087+ 0.001 


Ql 


1.084+ 0.044 


0.12+ 0.03 


71.4+ 14.6 


0.148+ 0.007 


Q2 


1.766+ 0.052 


0.20+ 0.01 


135.4+ 2.9 


0.098+ 0.001 


43 GHz -2001.20 


C 


3.733+0.149 






0.120+ 0.006 


Q2 


2.001+0.080 


0.19+ 0.03 


135.9+9.5 


0.114+ 0.006 


43 GHz -2001.29 


C 


3.097+ 0.104 






0.089+ 0.001 


Ql 


0.662+ 0.029 


0.15+ 0.03 


73.7+ 12.4 


0.142+ 0.007 


Q2 


1.815+0.061 


0.20+ 0.01 


135.8+ 2.9 


0.084+ 0.001 


43 GHz -2001.38 


C 


2.781+0.107 






0.089+ 0.001 


Ql 


0.439+ 0.021 


0.17+ 0.03 


66.3+ 10.7 


0.107+0.005 


Q2 


2.159+0.083 


0.19+0.01 


136.8+3.0 


0.089+ 0.001 


43 GHz -2001.48 


C 


2.692+ 0.104 






0.091+0.001 


Ql 


0.344+ 0.016 


0.18+0.03 


68.3+ 10.2 


0.065+0.003 


Q2 


2.219+ 0.086 


0.19+ 0.01 


138.9+3.0 


0.097+ 0.001 


43 GHz -2001.62 


C 


2.623+ 0.123 






0.082+ 0.001 


Ql 


0.456+ 0.025 


0.17+ 0.03 


67.6+ 10.6 


0.108+ 0.005 


Q2 


2.260+ 0.106 


0.20+ 0.01 


140.1+2.9 


0.096+ 0.001 






43 GHz -2001 


.70 




C 


2.440+ 0.159 






0.089+0.001 


Ql 


< 0.250 


0.19+0.03 


67.6+ 9.6 


<0.030 


Q2 


2.068+0.135 


0.20+0.01 


139.6+ 2.8 


0.093+0.001 


43 GHz -2001.78 


C 


3.036+0.173 






0.095+0.001 


Ql 


0.437+ 0.028 


0.20+ 0.03 


71.5+9.3 


0.106+ 0.005 


Q2 


2.184+ 0.125 


0.21+0.01 


137.7+ 2.7 


0.090+ 0.001 


43 GHz -2001.86 


C 


3.600+ 0.162 






0.054+ 0.001 


Ql 


0.594+ 0.032 


0.20+ 0.03 


84.4+ 9.2 


0.186+ 0.009 


Q2 


1.625+0.073 


0.25+0.01 


138.4+ 2.3 


0.063+ 0.001 


Q3 


0.944+ 0.042 


0.15+0.01 


117.6+3.9 


0.031+0.001 


43 GHz - 2002.07 


C 


4.334+ 0.145 






0.086+0.001 


Ql 


0.790+ 0.035 


0.16+ 0.03 


83.9+ 11.4 


0.201+0.010 


Q2 


1.508+0.050 


0.26+ 0.01 


143.0+ 2.2 


0.100+ 0.001 


Q3 


< 1.069 


0.19+0.01 


125.6+3.0 


<0.030 
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Table 2. continued. 



Comp. 


S 


r 





FWHM 




(Jy) 


(mas) 


O 


(mas) 


43 GHz - 2002.40 


C 


4.204+ 0.177 






0.046+ 0.001 


Qi 


0.884+ 0.045 


0.18+0.03 


92.9+9.9 


0.298+0.015 


Q2 


0.921+0.039 


0.28+ 0.01 


144.4+ 2.0 


0.032+ 0.001 


Q3 


1.695± 0.086 


0.17±0.03 


127.8+11.0 


0.102+ 0.005 


43 GHz - 2002.48 


C 


4.409± 0.171 






0.058+ 0.001 


Qi 


0.576+ 0.028 


0.21+0.03 


86.6+ 8.7 


0.216+ 0.011 


Q2 


1.045+0.040 


0.28+0.01 


144.0+2.1 


0.060+ 0.001 


Q3 


1.540+ 0.074 


0.17+0.03 


129.8+ 11.1 


0.103+0.005 


43 GHz -2003.20 


C 


3.356+0.145 






0.085+0.001 


QI 


0.582± 0.030 


0.23± 0.03 


101.0+8.1 


0.180+ 0.009 


Q2 


0.788± 0.034 


0.32± 0.01 


150.8+ 1.8 


0.050+ 0.001 


Q3 


1.440± 0.074 


0.19+0.03 


146.7+ 9.9 


0.107+ 0.005 


43 GHz -2003.49 


C 


4.647+ 0.259 






0.088+ 0.001 


QI 


0.554+ 0.035 


0.24+ 0.03 


95.6+7.7 


0.204+ 0.010 


Q2 


< 0.569 


0.31±0.01 


151.5+ 1.8 


<0.030 


Q3 


1.596± 0.089 


0.19± 0.01 


148.5+3.0 


0.097+ 0.001 


43 GHz - 2004.64 


C 


4.937+0.160 






0.107+0.005 


QI 


0.558+0.018 


0.30+ 0.03 


108.0+ 6.1 


0.229+0.011 


Q2 


0.438± 0.007 


0.38±0.01 


155.6+ 1.5 


0.035+ 0.000 


Q3 


1.063± 0.017 


0.23+0.01 


151.8+ 2.5 


0.098+ 0.001 


43 GHz - 2004.97 


C 


3.229+0.116 






0.128+ 0.006 


QI 


0.475+ 0.017 


0.31+0.03 


110.6+6.0 


0.204+ 0.010 


Q2 


0.507+0.018 


0.40+ 0.03 


158.1+4.6 


0.069+ 0.003 


Q3 


1.160+0.042 


0.21 + 0.03 


157.6+ 8.6 


0.117+0.006 


43 GHz -2005.25 


C 


2.512+ 0.093 






0.119+0.006 


QI 


0.435+0.016 


0.33+0.03 


111.7+5.6 


0.226+0.011 


Q2 


< 0.875 


0.40+ 0.01 


160.3+ 1.4 


<0.030 


Q3 


1.009± 0.037 


0.22+ 0.03 


162.3+ 8.6 


0.153+0.008 


43 GHz -2006.17 


C 


1.280± 0.049 






0.148+ 0.007 


QI 


0.261+0.010 


0.36+ 0.03 


115.6+5.1 


0.217+0.011 


Q2 


1.788+0.045 


0.36+ 0.01 


160.2+ 1.6 


0.039+ 0.000 


Q3 


0.940± 0.024 


0.25± 0.01 


164.8+ 2.3 


0.087+ 0.001 


43 GHz - 2006.24 


C 


1.257+0.049 






0.142+ 0.007 


QI 


0.349± 0.014 


0.34± 0.03 


113.7+5.3 


0.233+0.012 


Q2 


< 1.484 


0.36+ 0.01 


160.4+ 1.6 


<0.030 


Q3 


1.218+0.033 


0.26+ 0.01 


163.5+2.2 


0.082+ 0.001 


43 GHz - 2006.30 


C 


1.207+0.049 






0.153+0.008 


QI 


0.362± 0.015 


0.34± 0.03 


112.1+5.3 


0.248+ 0.012 


Q2 


1.470± 0.043 


0.35+0.01 


160.2+ 1.6 


0.051+0.001 


Q3 


1.295± 0.038 


0.25+0.01 


163.5+2.2 


0.077+ 0.001 


43 GHz - 2006.37 


C 


1.133+0.047 






0.155+0.008 


QI 


0.390+ 0.016 


0.33+0.03 


115.8+5.6 


0.262+ 0.013 


Q2 


1.236± 0.037 


0.35+0.01 


160.6+ 1.6 


0.056+ 0.001 
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Table 2. continued. 



Comp. 


S 


r 





FWHM 




(Jy) 


(mas) 


o 


(mas) 


Q3 


1.430± 0.043 


0.25± 0.01 


163.4± 2.3 


0.077± 0.001 


43 GHz - 2006.46 


C 


1.079+0.045 






0.149+0.007 


Ql 


0.338+0.014 


0.35+0.03 


114.6+5.3 


0.258+0.013 


Q2 


1.038+0.032 


0.36+0.01 


161.0+ 1.6 


0.049+ 0.000 


Q3 


1.742± 0.054 


0.25+0.01 


162.8± 2.3 


0.079± 0.001 


43 GHz - 2006.88 


C 


0.760+ 0.027 






0.150+ 0.008 


Ql 


0.253+ 0.009 


0.38+0.03 


117.5+4.8 


0.274+ 0.014 


Q2 


0.769± 0.017 


0.42+ 0.01 


163.2+ 1.4 


0.049+ 0.000 


Q3 


2.465± 0.055 


0.24± 0.01 


163.4± 2.4 


0.076± 0.001 


43 GHz -2007.16 


C 


0.903 ± 0.052 






0.170+ 0.009 


Ql 


0.348+ 0.020 


0.38+ 0.03 


119.7+ 4.8 


0.258+0.013 


Q2 


1.143+0.057 


0.45+0.01 


163.8+ 1.3 


0.085+0.001 


Q3 


3.895± 0.222 


0.23+ 0.03 


164.0+ 7.9 


0.101+0.005 


8.4 GHz - 2002.48 


Core 


6.352± 0.635 






0.378± 0.038 


X4 


i 0.088 


0.62+ 0.24 


-98.8+22.3 


<0.180 


X5 


0.073+ 0.007 


1.73+0.42 


-52.5+ 13.9 


0.836+0.084 


X6 


0.116+0.012 


3.59+ 0.44 


-43.0+7.1 


0.883+0.088 


X7 


0.103+0.010 


24.71+ 1.39 


-31.8+3.2 


2.785+ 0.278 
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Table 3. Proper motion results. For each model component, we give the number of data points for the fit of the trajectory (N), the mean position 
angle (< >), the mean angular speed (< >), the mean proper motion (< ju >, and its corresponding jS^pp< ju >), as well as its radial (< jUrad >) 
and normal component (< /x„orm >. and its corresponding jS^p_norm)- 



Comp. 


N 


<0> 


<&> 


<JU> 


oobs 
A'app 


<J«rad > 




oohs 

"app_norm 








7yr 


mas/yr 


C 


mas/yr 


mas/yr 


C 


Ql 


33 


139.5+7.1 


10.7±0.7 


0.047+0.002 


3.26+0.14 


0.026+0.002 


0.039±0.003 


2.71+0.21 


Q2 


34 


175.4±1.3 


5.8+0.3 


0.042±0.001 


2.85+0.07 


0.036±0.001 


0.020±0.003 


1.39+0.21 


Q3 


16 


153.8+2.7 


8.0±0.5 


0.034±0.002 


2.29±0.14 


0.018±0.002 


0.029±0.002 


2.01±0.14 



